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Abstract: In the framework of the dipole picture of the BFKL pomeron we discuss a few 
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1. Introduction 



Testing the BFKL pomeron through collisions of tagged e + e~ pairs with very large mo- 
mentum transfers is an attractive possibility which has already been discussed in Refs. [|]. Q . 
Clearly, the crucial ingredient in the expression for the e + e~ total cross section at fixed 
Q\, Q\ momentum transfers of the tagged leptons is the total cross section of two virtual 
photons of "masses" Q\ and Q\. Calculation of this cross section, <t 77 , is basic for the 
content of Refs. |, ||. 

Here, we present a method of calculating a 77 , alternative to the one presented in U 
and 0, which is an implementation of the dipole picture of the BFKL pomeron proposed 
in Refs. || f|, || This implementation goes somewhat beyond the applications of the 
dipole picture given in Refs. J?], |8|, |, [T0[ [□], [L2| and is based on a discussion of the role of 
highly asymmetric q — q configurations of light quarks in the onium - onium scattering at 
high energy within the framework of Mueller's QCD dipole picture fl3| . 



We start with the forward onium - onium amplitude for a single pomeron exchange, 
which we calculated ab initio. We believe that this new expression is better than the one 
used in Refs. || [LI], Its detailed derivation is given in Appendix A and the result is 



F« = rca 2 r A r B [ ^L e ^ Y ( ^V"* h(j) 
J 2ni KruJ 



r B J 

Here a is the strong coupling constant, N is the number of colors, and r B are the transverse 
sizes of the two colliding onia, A (7) = aNx{^)/^ where 



x(7) = 2^(i)-^(i-h)-^(h) 

and ^ 

7 2 (2-7) 2 

The quantity 

y = io g (-) (4) 

is the total length of the dipole cascade, i.e., the sum of the cascade lengths of the two 
colliding onia. s is the total cm. energy of the collision and Sq is the relevant scale of 
the problem. s cannot be calculated within the leading logarithmic approximation and 



if, 



dlogT 



(2) 
(3) 
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therefore it remains an unknown element in this approach. Its determination must rely on 
one's physical intuition and on results of a phenomenological analysis of data. In the present 
paper we explore the consequences of the choice suggested by the dipole picture 

While (HD reproduces the saddle point approximation of derived in [||, |J [5| (and 
employed in |9|, [TO], [TTj, 0), it also contains contributions which are neglected in the contour 
integral representations of F^ employed in |10 , [11], 12| . In other words, in Refs. |9|, 



PH , [l2| , these components of the integrand in the contour integral representation of F^ 



which become unity at the saddle point, are kept equal one throughout the whole contour 
integration. This approximation has been corrected in our present expression for F^\ 

The way to employ the dipole picture to calculate the total 7*7* cross section is, in 
principle, straightforward. From F^ 1 ' and the well known (compare, e.g., || [XOj, O, |14f ) 
wave functions of the two photons, A and B, of the virtual masses Qa,b, longitudinally (L) 
or transversely (T) polarized, {ta,b, za,b] Qa,b), we obtain the forward 7*7* amplitude 

F 77 = J \^ T (r A ,z A -,Q A )\ 2 \^ L ' T (r B ,z B -,Q B )\ 2 F^d 2 r A dzAd 2 r B dz B , (5) 
and the total cross section which, with our conventions, reads 

cr 77 = 2 Re F 77 . (6) 

Then, in order to evaluate the integral in (||), we have to decide what to take for Y, the 
length of the cascade. According to (|j) this amounts to a selection of the scale so- Two 
choices were discussed in the literature. In Ref. |l[ sq was taken as 

so = cQaQb , (7) 

with c = 100. This apparently natural choice has an attractive feature to be a simple analytic 
function of Qa and Q B . Another possibility was to take 

s = cQl, (8) 
where Q > is the larger of Qa , Qb- This gives 

Y = ]og(—) , (9) 
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a formula which provided an excellent fit to the proton structure function at small xbj 0, HQ- 



Following ||T3| we observe, however, that from the point of view of the dipole picture 
neither (^) nor (|j) is really satisfactory. The point is that Eq. (|U) refers to collisions of 
two dipoles and thus the relevant scale so must be expressed in term of the parameters 
characterizing these dipoles (i.e., longitudinal momenta za, zb and transverse sizes r^rg) 
rather than Qa and Qb- This is clear if one observes that Qa and Qb are not even defined in 
(|l|). The possible choices of s consistent with the dipole picture were discussed in |13| , where 



also a definite formula for Y was suggested. In the present paper we explore consequences 
of this choice for <r 77 and compare it with the results following from (|7]) and (H). We hope 
that our results shall be useful in testing the validity of the dipole picture approach in the 
small xbj physics. 

In the next Section we present our formulae for a 77 following from Y worked out in Ref. 
13| . In Section 3 we give our numerical results and their discussion. Section 4 contains the 



conclusions. Appendix A presents a detailed derivation of the formula (|IJ), in Appendix B 
we calculate <r 77 for three factorizable forms of Y defined by: (|T3"D, (0), and (|8|). The case (|7|) 
gives cr 77 of Ref. |J whereas the last one, (H), results in cr 77 which follows from the version 
of the dipole picture implemented in Ref. 

2. The total 7* — 7* cross section 

In this Section we derive the formulae for the total cross section of two virtual gammas 
using the formula ([[]) for onium - onium forward amplitude derived in Appendix A, with Y 



taken from Ref. [13 1 



' < < 1 2 

sz A z B r A' B 



Y = y A + yB = \og[ - A Z A ^ , (10) 



\ CT int J 

where c is the arbitrary constant of the leading log approximation, s = AEaEb is the square 
of the total cm. energy of the colliding virtual photons, 

, \ z if z < i , 

z < ={ ~ 2 11 

and 

r int = const r> (12) 



2 Note that @ and (g) give rather different results for Q 2 dependence of 7 — 7 cross section which should 
not be too difficult to test once the relevant data are available. 
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where r> is the larger of r A and tb- Tint is interpreted (see [Q) as the time needed for the 
exchanged gluons to travel the necessary distance in the transverse space. 

We will confront the results obtained with (|l~2"D with the ones one gets replacing ([12]) by 
a symmetric expression 

T L -» T 2 = r A r B (13) 

which leads to formulae close to the ones advocated in Refs. [1], 0] . 

To get 7* — 7* amplitude we employ now the wave functions of the virtual photons, 
^(ta, %a\ Qa), ^{tb, zb] Qb), and calculate F 77 of (f|), where for the transverse (T) and 



longitudinal (L) photons we have (compare || |10[ 0, |14 



Na em e 2 f 



\* T ' L (r,z;Q)\ 2 = $ T ' L (r } z;Q) = ^f±W T > L (r, z-Q) , (14) 



7T Z 



W T (r,z;Q) = \[z 2 + (1 - zf}Q 2 Kl{Qr) , (15) 



W L (r, z; Q) = 2z(l - z)Q 2 K 2 Q (Qr) , (16) 

where Q = J z{l — z)Q, a em = 1/137 and e"j = 2/3 (the sum of the squares of the charges 
of three quarks). 

Inserting fllPl) into (jl]) and employing (|5]) we obtain for the total 7*7* cross section 

1 

cx 77 = 4(2vr) 3 a 2 / dr A r\ dz A $ T > L (z A ,r A ;Q A ) 
Jo Jo 

x r dr B rl [~ 2 dz B ^ L (z B , r B ; Qb) [ ^.e^ Y (^) ^(7) , (17) 
Jo Jo J 2m \rs/ 

where Y = log£ and £ = {sz A z^r A r B /(cT 2 nt )). Note that since $ T and $ L are invariant 
against the replacements: z A)B — > (1 — z AjB ) and (1 — Za,b) ~^ z a,b, we can drop the < 

1 1 /2 

superscripts and integrate over z s as follows J dz — > 2 J dz . 

The expression (|T7| ) is comparatively easy to evaluate when, as in (|T3|), r^ t — > r 2 = r^re 
(the case close to the one of Refs. |], U) because the integrals factorize into integrals over 
the A and B variables, the integrations over r A and tb can be done analytically and we 
obtain 
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n QaQb-J 2iri\cQ A Q B J \Qa- 

where H L ' T = 4 A ^ Z l > t (-y)S l > t (^)Z l ' t (2 - 7 )5 L ' T (2 - 7) is the product of the functions of 
7 defined in the Appendix B. 

In the case of (0), however, there is no factorization and one has to face a 5-dimensional 
integration with one of the integrals being a contour integral in the complex plane 7 (along 
a straight line parallel to the imaginary axis). This forces us to use a numerical method for 
the evaluation of Eq. fllTD - It turns out (see the discussion of our results below) that in the 
very high energy limit (s/c very large) one can safely use the saddle point approximation for 
the contour integral 



^l e A( 7) y^y-\ (7) = lJ^ h{j0 )^ Pe -W^/^ , (19) 
2m KtbJ V ^ 

with Y = log£, £ — s/sq, with s as the case might be (see above), = [7aN((3) log(£)/7r] _1 

and A p = — 7o = 1 — a^log^/re) is the saddle point which, in the limit s/c — > 00, 

equals 1. Then the numerical integration reduces to 4 dimensions. Note that ([19]) exhibits 

the source of the substantial difference in dependences on Qb/Qa following from (|7|) and 

(181), see Fig. 2. This difference sits in the £ Ap factor: 



UasJ agamst WJ ' (20) 

3. Numerical results and discussion 

We considered 4 cases of Y = log(s/s ) and calculated the corresponding cross sections: 
(a) The case of Eq. (0), s = c^)QaQb, employed in Refs. Jl], U, 



(b) the case of Eq. (||), s = C(b)Q>, employed in Ref. [[/], 

(c) the case of Eq. flT2]), s = (c^rl > /(z A z^r\r^)), discussed in Ref. [ITS], 

(d) the case of Eq. fljjD , s = (c^d) I (zA z B r A r b)) , discussed also in Ref. fT3[ . 



Comments: As shown in Appendix B, in the Case (a) we obtain the same formula for <r 77 
as in Ref. [0. Also let us note that, when Qa = Qb and the arbitrary constants are set to 
the same value C( a ) = c^) = c, the cr's for Case (a) and Case (b) are identical. 
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The Cases (a) - - (d) were calculated in the saddle point approximation given by the 
formula fll9|) and subsequent 4-dimensional integration. In Cases (a), (b) and (d), we checked 
the accuracy of this procedure calculating cx 77 analytically up to the final contour integration 
over 7 which was done with the help of MATHEMATICA. It turned out that the results of 
these two procedures agree to within 15 percent. 

In order to exhibit the asymmetries when Qa ^ Qb we introduced the asymmetry pa- 
rameter, (, defined as 

C = ^. (21) 

From the asymptotic forms ( [56] ) and ([54]) we see that the asymmetry in Qa, in Cases (a) 
and (d) is given approximately by the factor 

e -fa e log 2 (C)_ ( 22 ) 

In the asymmetric cases, especially in Case (b), this estimate is not good enough. 

Clearly, the choices of the values of the arbitrary constant c involved in all Y's discussed 
in this paper are very important in determining the size of the cross section. They can either 
be fitted to experimental results (compare Ref. [^]) or set following some prejudices of the 
authors (compare, e.g., Ref. M): in Ref. [I| 



and in Ref. 00 



c = C(a) = 100, £ = £(a) = (23) 

c = c (b) = 0.57, £ = = — 8 — . (24) 

The constants c for Cases (c) and (d) were set to fit the ex's of Cases (b) and (a), respectively, 
for Q A = Q B = 4 GeV and y/s = 200 GeV. They come out to be: C( c ) = 0.0055, and C(d) = 
2.5. 

To estimate the role of c's it is enough to use the asymptotic formula ( |64]) for Qa = Qb- 
We get 



~K\ = it) tP" < 25 > 
^77(0) U(6)/ V%» 



3 The authors are grateful to R. Peschanski for providing them with C(&) of (]2J) which gives the fit of Ref. 
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where is given below (|19|). Since, in the limit s — > oo, 0(log£( o )) = (log £(&)), we have 
approximately 



= fMf p . (26) 

cr 77 (o) Vc (a )/ 

In Fig. 1 we present the cr 77 's for Cases (a) - - (d), for Qa = Qb, setting the strong 
coupling constant a = 0.11, hence A p = 0.3. The values of the constants c were taken as 
in Ref. Q (c( Q ) = 100) and @ (c(&) = 0.57). The resulting cross sections differ appreciably, 
consistently with Eq. (|26|). The dipole model results, (c) and (d), were fitted to the predic- 
tions of (b) and (a), respectively, at the point Q 2 = 16 GeV 2 and y/s = 200 GeV. One sees 
that they follow closely the results of (a) and (b) for all considered values of Q 2 and y/s. 




Figure 1: Total cross section in the dominant TT channel for Qa = Qb = Q- The solid lines 
represent the results for Cases (b) and (c), whereas the dashed lines show Cases (a) and (d). Using 
the values of the scale parameter c given in the text, one finds that Cases (b) and (c) coincide in 
the wide range of Q 2 . The results for Cases (a) and (d) slightly differ for very large Q 2 . 
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In Fig. 2 the dependence on the asymmetry parameter ( = Qb/Qa is plotted. Other 
parameters are chosen as in Fig. 1. One sees that the ^-dependence is almost identical for 
two versions of the dipole model ((c) and (d)) and the symmetric proposal of 0, @j- The 
case (b) differs significantly, however, from the others, giving a much stronger dependence 
on C 




Figure 2: Total cross section in the dominant TT channel for 1 < £ < 6. 
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4. Conclusions 



The predictions of the dipole model for the photon - photon cross section depend strongly 
on the scale determining the length of the dipole cascade in the incident photons. The scale 
suggested previously || gives substantially smaller cross section than the one suggested 
by a fit of the dipole model results to the proton structure function J7[. On the other hand, 
the dependence of the cross section on the ratio Qb / Qa fo r the two colliding photons turned 
out to be the same for the two extreme cases of the dipole model ((c), (d)) suggesting that 
it hardly depends on the details of the model. 

We conclude that future measurements of 7*7* cross section may be useful in determining 
the length of the dipole (gluon) cascade but, probably, not very helpful in understanding the 
details of the dipole-dipole interaction. This makes rather urgent the need of determining 
the relevant scales from the higher-order perturbative calculations. 

Acknowledgment: The authors thank W. Broniowski for his helpful suggestions concerning 
the usage of the Monte-Carlo method. 

Appendix A. 

Forward onium - onium amplitude including the light-cone momentum fractions 
dependencies 

In this Appendix we give a direct derivation of the formula for = J d 2 bF^ l \b) in the 
form of a contour integral. This derivation accommodates the option that the dipole - dipole 
cross section takes the asymptotic form only above some energy Eq of the slowest gluon. 

We follow the notation of the main text: the subscripts A and B refer to one of the 
virtual photons. For example, xa is the transverse size of a dipole originating from the 
virtual photon A. Also x < denotes the smaller of the two sizes xa and xb- As before, r A ,B 
and za,b are the arguments of the photon wave functions. Remember, however, that z A B 
are defined by (|TT|). Consulting of Ref. [[OJ may also be helpful. 

We start from the formula (3) of Mueller's paper || which we write in the form 

/dec doc q 
d 2 sa(x A , x B )n(r A , x A , y A , s)n(r B , x B , Vb, b-s) , (27) 
x A x B 

where in the limit of high energy 

r dl 

a{x A ,x B ) =4 J -[1- J (x A l)][l- Mx B l)\ =4[l + log(x>/x<)] , (28) 
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and 

n(r } x } y } s) = ^J Vfl-^f^V.Vj) , (29) 

with 



I(r,x, 5,7) = / d 2 «; + \x — w\\s — \x — w\Y (|~r — w\\ — ~r — u>|) 7 . (30) 
Following Mueller we also find (see 



z < 



2/ = log — , (31) 

V Zq J 

where z K is defined in ( |TTD and Zq is the minimal fraction of the onium energy which can be 
carried by a dipole in the wave function (actually the minimal energy carried by one of the 
gluons forming a dipole). Mueller does not explain why the integration must stop at some zq 
(except for the obvious fact that for z = the formula diverges). It is clear, however, that 
without understanding this problem there is no chance to estimate reasonably zq. Following 
the arguments of |T3j {e.g., Eqs. (18) and (21)) we take 

_ const r> _ p£ 

zo — — + ~~ ~r ' y 6Z > 

where p$ = 2E , E Q being the energy of the slowest gluon. 
When all this is substituted into ([!]) we obtain 



(1) a 2 fd lA d lB / 4g|rj r mj n y( z»p%rl \^B) 

= ^J 2™ 2™ {-^-- ) (1 - 1A) {-^—) (1 " 1B) ^ ^ ^ 1B) ' 

(33) 



where 



tt{r A , Tb, 1a, 1b) = [ X A XB ct{xa, x b ) ( — ) (— ) I(x A , x B , r A , r B , 1a, 1b) , (34) 

J X A X B \X A J \X B J 

with 

I(x A ,x B ,r a ,r B ,lA,lB) = J d 2 bd 2 sd 2 wd 2 u 

Mo — s + ^x B — u\\b — s — ^x B — u\\ \\2 r B ~ u \\~ 2 Tb ~ u \) 
\s + \x A — w\\s — \x A — w\y A (\^r A — w\\ — ^r A — w\\ 1A . (35) 



X 
X 



After the change of variables, s — w —> s, b — s — u — > b, Q35|) factorizes: 
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I(x A , x B , r A , r B , 7 A , nf B ) = I(x B , 2 - i B )I{x A , 2 - i A )I{r A , i A )I{r B , j B ) 



where 



/(x, A) = I d 2 s{\s + \x\\s — \x 



2' 



/(x, A) can be calculated using the technique of Mueller. Here we give only the 
reads 



where 



Consequently, we can write 



H{\) 



I(X,\) = 7TX 2{1 " X) H(\) , 

r 2 (i-iA)r(A-i; 



r 2 (|A)r(2 - A) 



V(r A , r B , 1a, 7b) = ^H(2 - lA )H(2 - 7B )#(7 A )#( 7B )r A - 7 V|- 7 M7A, 7b 



where 



^(1a,1b) 



00 dx A dx B 7A _ 2 7B _ 2 ( v 



/o X A X B 

Integration over x A and x B is tedious but we do it explicitly 

^ dx 



w(7a,7b) 



o Xa 



X A 



(x B ) 7s (l + log — ) 

X B X B 



+ x A r^(x B y*- 2 (i + iog Xi; 



'XA %B 

The integrals over dx B are well defined and are: 

CA dx B . . „ 



X A 



r 2 
Xa 



XB 

00 dx B . 



IbT 
3 - 7b 
(2 - 7b) 2 

When this is substituted into (03) and integrated over c/xa we obtain 



x A x B 



k xb) 



1 + log 


za\ 


— ™7B 1 

— X A ^ 




x B ) 


'l + lo { 


X B ^ 

T 


\ _ ™7B 
J - X A ( 







^(7a,7b) = h( 7 



Bj 



dx 



XA 



Xa 



Ki B )t 



7.4+7S-2 



2 - 7a - 7b 
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where we have regularized the integral by taking the lower limit to be e. h{pf B ) is given by 



Inserting (|39| ) and ( f4"5D into P0| ) we obtain 

(l-7A) 2 (l-7 B mrA,r B , 7 A,7 B ) = vrV^r^^^ . (47) 

2 - 7a - 1b 

Now, it turns out that the result for FW depends on the order of integration over 7's. The 
integral over j B in (|33|) can be calculated just by taking the residue at the pole (2—7^4— 7_b)~ 1 . 
We obtain 

F« = ito?r A r B f ^e^ Y (^) ^h(j) , (48) 
J 2m \r B J 

where 



Y = logip+p+z^rl/c) = log(s44 rj/c) , (49) 

and s = p^Ps i s the total cm. energy squared. The integration over 7^ can be similarly 
executed and one obtains for again the result ([48]) but with r A ,r B — > r^,^ exchanged. 
However, (|48|) is in fact invariant under this exchange (see, e.g., the asymptotic form of F^ 
(|T9"|) which exhibits such an invariance). 



Appendix B 



As it was pointed out in jTjJ, the dipole picture of BFKL ||, f|, Q gives some freedom 
of choice of the specific realizations of Y in ([[]). In this Appendix we give derivations of 
<j 77 for two cases of the explicitly symmetric (with respect to the virtual photons A and B) 
expressions for Y. 

We start with the case of Tj? t given by Eq. ([13]). Now F^ becomes 

FW = 7ia 2 r A r B J ^( S r A r B z<z</c)^{^y' l h{ 1 ) , (50) 

From (H), (|i~4]), (^) and (fl6|) it follows that the integrations over r A , r B , z A , and z B factorize. 
The integrals over r A and r B can be done with the help of the formula 6.576/4 of Gradstein 
and Ryzhik [15 ], and one obtains after a few straightforward operations 
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where 



and 



x ^ T ( 7 )^ T ( 7 )Z^(2- 7 )^ T (2- 7 ) , (51) 



S T h) = ^S L (7) , S%) = |7 ; } , (52) 



2-7- W/ ' W/ r(4- 7 _) ' 



Z T ( 7 ) = ( 2 dz\z 2 + (1 - z) 2 ]zi 7+-1 (l - z)h-- x , (53) 

JO 



l 

Z L ( 7 ) = 4 / 2 d^ 7+ (l - z)h~ , (54) 
Jo 



7± = 7 ± A( 7 ) . (55) 

It is also of interest to have the asymptotic (i.e., s/c — > oo) expression for cx 77 . This 
can be done through the saddle point approximation around 7 = 1 (from (Q) we have that 
X'(l) = 0). We obtain 



where 



a 77 = -{aNa em e)Y-^ x \^^e-^^ , (56) 

7T QaQb V 7T 



Ap = A( 7 = l), a c =[7aJVC(3)log(0/7r]- 1 , (57) 



cQ A Q 



B 



s T (i) = s L (i)^^ , g^(i) = V2 2 ; y , (5^ 

v ; v J l + A P ' v ; T(3 + A P ) ' v 



Z T (1) = [ 2 dz[z 2 + (1 - z^zh^-D^ _ z )-i(Ap+i) (59) 

JO 

1 

Z L (1) = 4 / 2 c^ (1+Ap) (l - ^ (1 ~ Ap) . (60) 
Jo 
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To have a direct comparison of the dipole picture formulae with those of Refs. |], ^[ and 
with Ref. 0, we employ the forms of Y which do not depend on r's and z's (((|) and (0), 
(U) and (H)). Now all integrals over ta , re, za, Zb can be done analytically employing the 
formula 6.576/4 of Ref. fl5| for the integrals over r's, and the integrals over z's are now 



simple Euler f3 functions. We obtain 
where 

l7 ' 2-7 r(2 + 7 ) r(4- 7 ) ' 1 ; 

^ ™ " 2 r(2 + 7) r(4 - 7) ' 11 

Setting £ = s/(cQaQb) we obtain <r 77 of Ref. [l|]. The transition between Eq. (8) of Ref. 
|TJ and Q5"l"| ) goes through the substitution 7 = 27' which changes \ of (0) and transforms 
the r.h.s. of ( |BTD into the r.h.s. of Eq. (8) in multiplied by a factor |. This factor can 
be traced to an approximation made in Ref. || . 

For the sake of completeness we give also the saddle point approximation of ( |61~D 

<t 77 = , (64) 

ttQaQb 




where 

r(l) = - ,1^(1) = — . (65) 
v ; 256 ' V ' 256 K J 

Note that for the case discussed in Ref. J7| we have £ = s/(cQ?J (compare Eq. (H), and the 

same set of formulae for <r 77 (|6~1~D - (164]) ). 



15 



References 



[1] S.J. Brodsky, F. Hautmann and D.E. Soper: Phys. Rev. Lett. 78 (1997) 803, see also, 
S.J. Brodsky, in Proceedings of the Workshop on Physics with High Energy Colliders, 
Brookhaven National Laboratory, Upton, New York, 1996, edited by Peter Jamin (to 
be published), and F. Hautmann, Proceedings of the XXVIII International Conference 
on High Energy Physics, Warsaw, 1996, edited by A. Wroblewski (to be published) 

[2] J. Bartels, A. De Broek and H. Lotter: Phys. Lett. B389 (1996) 742 

[3] A.H. Mueller: Nucl. Phys. B415 (1994) 373 

[4] A.H. Mueller and B. Patel: Nucl. Phys. B425 (1994) 471 

[5] A.H. Mueller: Nucl. Phys. B437 (1995) 107 

[6] N.N. Nikolaev and B.G. Zakharov: Z. Phys C49 (1991) 607; C53 (1992) 331; C64 
(1994) 651; J. Exp. Theor. Phys. 78 (1994) 598 

[7] H. Navalet, R. Peschanski and Ch. Royon: Phys. Lett. B366, (1995) 329 

[8] H. Navalet, R. Peschanski, Ch. Royon and S. Wallon: Phys. Lett. B385 (1996) 375 

[9] A. Bialas and R. Peschanski: Phys. Lett. B378 (1996) 302 

[10] A. Bialas and R. Peschanski: Phys. Lett. B387 (1996) 405 

[11] A. Bialas: Acta Phys. Pol. B27 (1996) 1263 

[12] A. Bialas, W. Czyz and W. Florkowski: Phys. Rev. D55 (1997) 6830 

[13] A. Bialas: Acta Phys. Pol. B28 (1997) 1239 

[14] J. Bjorken, J. Kogut, D. Soper: Phys. Rev. D3 (1971) 1382 

[15] P.S. Gradstein and I.M. Ryzhik: Tablicy integralov, sum, ryadov i proizvedeni Moskva, 
1971, 5th edition, p.707, formula 6.576/4 



